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Abstract 

The hysteresis properties of ferromagnetic materials at low field are described by 
the Rayleigh law. We analyze the problem in the light of modern statistical mechan- 
ics models of hysteresis. In particular, we compute the demagnetization curve and 
derive the Rayleigh parameters a and h in the random-field Ising model and in a 
model of domain wall depinning. In the random-field Ising model the Rayleigh law 
is obeyed only in the disorder dominated phase, while in the low disorder phase it is 
not possible to demagnetize the sample. This approach allows us to link a and b to 
microstructural parameters, such as the domain wall energy, the internal disorder 
or the exchange interactions. Finally, our results are compared with experiments. 
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1 Introduction 



In 1887 Lord Rayleigh analyzed the hysteresis properties of ferromagnetic 
materials at low fields, close to the demagnetized state [1-3]. When the field is 
cycled between —Hm and , the magnetization M is found to follow a simple 
quadratic law M = (a + bHm) ± b{H^ — H'^)/2, where the signs ± distinguish 
the upper and lower branch of the loop. The area of the loop can easily be 
computed and is given hj W = (4/3)6if^. In addition, the response of the 
system to a small field change, starting from the demagnetized state is given 
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by Mm — ai?m =t bH^. The Rayleigh law has been shown to hold in several 
ferromagnetic materials [4] and it has been also widely applied to describe 
ferroelectric materials [5,6]. A few papers have reported significant deviations 
from the simple quadratic law but no explanation has been provided [7] . 

In 1942 the Rayleigh law was interpreted by Ncel in terms of the motion 
of a point (i.e. a rigid domain wall) in a random energy landscape, whose 
statistical properties determine the value of a and b [8]. In particular, a is 
associated to reversible motions inside one of the many minima of the random 
potential, while b describes irreversible jumps between different valleys. Several 
generalizations of this approach have been proposed, considering more refined 
forms for the random energy landscape [9-13]. These approaches assume a 
single rigid domain wall and thus do not consider the effect of domain wall 
bowing and the interactions between different domains. When this effects are 
important it is hard to link the Rayleigh law to the microstructural properties 
of the material, such as dimensionality, chemical composition or interaction 
types. 

Here, we reconsider all these problems in light of novel statistical mechanics 
approaches to hysteresis. In particular we compute the Rayleigh parameters 
a and b in the random-field Ising model [14] and in a model of domain wall 
depinning [15-21]. This approach allows us to link a and b to microstruc- 
tural parameters, such as the domain wall energy, the internal disorder or the 
exchange interactions. 



2 Models 

A microscopic theory of hysteresis should in principle recover the Rayleigh law 
from the collective properties of interacting magnetic moments. The Rayleigh 
parameters would then be expressed in terms of micromagnetic parameters, 
that are available in the literature for a variety of ferromagnetic materials. 
Ferromagnets can in general be described by a locally varying magnetization 
Miif) evolving under the action of the externally applied field. The evolution 
of Mi is ruled by appropriate relaxation equations that can be written in 
terms of the micromagnetic energy [2,3]. Considering for simplicity a uniaxial 
material, we can write the energy as 

^ = E / dh[A{VMif + K{Minif - fioMsiH^ + H^^JM^ (1) 

where A is the exchange interaction, K is the anisotropy constant, is the 
anisotropy axis, and H^l^ is the component i of the demagnetizing field and 
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Ms is the saturation magnetization. 

An essential contribution to the properties of the hysteresis loops is given by 
disorder, due to crystal imperfections, internal stresses, non magnetic impuri- 
ties, that are present in most magnetic materials. Neel recognized this fact and 
replaced the micromagnetic free energy by a random function of the magneti- 
zation (i.e. E = Er{M) — HM) [8]. The approximation is rather drastic but 
makes the problem analytically tractable. In general we can model the disorder 
by quenched local fluctuations of the micromagnetic parameters. Disordered 
interaction terms are conventionally denoted as random bonds {A — > A{r)), 
random anisotropics (rij — > ni{r)) and random fields {H Hext + h{r)). 
A complete solution of micromagnetic equations including disorder is a very 
complicated task. One should then resort to some kind of approximation. In 
recent years two main approaches have been undertaken to describe the mag- 
netization properties of disordered materials. 

The first approach takes into account the fact that in soft magnetic materials 
demagnetizing fields give rise to broad domains and the magnetization process 
is dominated by domain wall motion. One can thus reduce the problem to the 
motion of a flexible domain wall in a random potential. The domain wall 

contribution to the micromagnetic free energy can be expressed in term of 
the domain wall coordinates z{x), inserting in Eq. (1) Mi{f) = 5i^Msg{{r^ — 

z{x))/5u,) in Eq. (1), where 5^ ~ \j AjK is the domain wall width and g{x) ~ 
±1 for a; ^ ±oo [15,18,21]: 

= / d^xYi^{Vz{x)f - zix)iioMs{H + Hai{z{x)}) + V{x, z{x))] (2) 

where the domain wall energy is given by 7^ ~ V AK, the stray field is given 
hj Hd = fioMg J rf^x' (?2;/(?a;'(a;—a;')/|r—r'|^ and V is a random function taking 
into account all the disorder contributions. In the following we will consider 
y as a superposition of pinning centers randomly distributed in space 

V{x,z) = ^/oexp-((^ - Z,{x))/Q' (3) 
p 

where /o is the strength of each pinning center, of width C,p located in Zp. The 
effective form of the pinning potential is not essential as long as the interaction 
range is finite. This is not the case for pinning due to isolated dislocations, 
whose stress field decays as 1/r. One should notice, however, that dislocations 
are typically arranged in patterns where the stress field is screened, providing 
us with an effective correlation length. In the following, we will restrict our 
attention to cases in which the pinning correlation length is smaller than the 
domain wall width, so that we can replace by 5w This case corresponds 
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mainly to the effect of non-magnetic impurities (see Ref. [16] for a discussion 
of this point) . 

Models based on Eq. 2 have been used in the past to compute the coercive 
field [15,17], analyze thermal relaxation [16], explain the statistical properties 
of the Barkhausen noise [19-22] and describe magnetization creep in thin films 
[23]. 

An approach based on domain wall motion is not adequate to describe the 
magnetization properties of hard ferromagnets, where the presence of strong 
random anisotropies prevents the formation of extended domains. For this 
class of materials a description in terms of interacting spins seems more ap- 
propriate. Several disordered spin models have been proposed in the past to 
describe hysteresis. Among those, the simplest and most studied is the RFIM 
[14,24], where a spin Sj = ±1 is assigned to each site i of a rf— dimensional 
lattice. The spins are coupled to their nearest-neighbors spins by a ferromag- 
netic interaction of strength J and to the external field H. In addition, to 
each site of the lattice it is associated a random field hi taken from a Gaussian 
probability distribution with variance R, P{h) — exp(— /i^/2i?^)/-\/27rit!. The 
Hamiltonian thus reads 

E^-Y^Js,s^-Y.{H + hi)s,, (4) 

where the first sum is restricted to nearest-neighbors pairs. The dynamics 
proposed in Ref. [14] is such that the spins align with the local field 

Si = sign(J^Sj + (5) 

3 

In this way a single spin flip can lead the neighboring spins to flip, eventually 
trigger an avalanche. Using this dynamics it has been shown that the RFIM 
displays a phase transition in d > 3 as a function of R [14,24]. For R < Rc, 

the saturation loop has a discrete jump in the magnetization at H = ±Hc. 
The jump disappears for R > Rc and R = Rc corresponds to a critical point 
and the model satisfies scaling laws [14,24]. 



3 Collective pinning effects in domain wall hysteresis 

The hysteresis properties of interfaces in random media have been studied 
in the past in the context of frictional sliding [25]. Similar studies for ferro- 
magnetic domain walls has only been restricted to the case of high driving 
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frequencies and large fields [26]. It is interesting to note, however, that in- 
terfaces in random media obey return point memory, a typical properties of 
ferromagnetic hysteresis. Here, we use a model based on Eq. (2) to analyze 
the Rayleigh law at low frequency. For the sake of simplicity we do not con- 
sider here demagnetizing fields, which are essential, however, to account for 
the large scale behavior of the magnetization [21]. When considering small 
scale displacements of the domain wall, we expect the domain wall energy 
(scaling as in Fourier space a deformation of wave vector q.) to be more 
relevant of the stray field contribution (scaling as q). The following arguments 
can nevertheless be generahzed including the effects of demagnetizing fields 
[15]. 

The scaling properties of the Rayleigh parameters can be obtained using col- 
lective pinning theory [15,16,23]. The central concept is the identification of a 
coherence length Lg, defining a region of space where the domain wall moves 
freely from the pinning centers. The coherence length can be obtained com- 
paring the domain wall energy with the disorder fluctuations over a region of 
length Lc, considering only small transverse wall deformations of the order of 

The depinning field He, which can be identified with the coercive field, is 
obtained comparing the pinning energy with the magnetostatic energy over a 
region of length [15,16,23] 

IIoM,H,LX ~ MnoLXy/' H, = (no/J)/(7«,5^/ioM,). (7) 

This expression for the coercive field recovers the result obtained in Ref. [15] 
when two dimensional domain vaulting is considered. In models based on a 
rigid domain wall [8,10,11] the coercive field is instead proportional to the 
standard deviation of the pinning field (i.e. He oc fo-y/no). 

Under the application of a small external field the "unpinned" regions of the 
domain wall will bow slightly and contribute to the magnetization [16]. The 
expression for the susceptibility a = dM/dH is similar to the one found in 
textbooks when discussing reversible susceptibility [2,3], with exception that 
the bowing length is now given by Lc, 

a ~ pioM^Ll/^^ ~ {lwSl)/{nof^). (8) 

This result differs from estimates based on generalizations of Neel theory 
[8,10,11] where ^yu, does not appear. 
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The first correction to the hnear susceptibihty, and hence the parameter b, 
can be obtained noticing that a smaU apphed field can in principle lead to 
local depinning events, thus increasing the coherence length [25]. Using the 
arguments reported in Ref. [25], one can show that, to lowest order in H, 
Lc{H) — Lc(0)(l + c\H\/Hc) where c is a numerical factor. Inserting this 
expression in Eq. (8) and expanding for small H, we obtain b ~ a/ He- A 
similar expression was already reported by Neel in 1942 [8]. 

In order to test these consideration, we perform numerical simulations based on 
Eq. 2. The coordinates of the domain walls evolve according to an overdamped 
equation Tdz/dt = —dE/6z, where T is an effective viscosity. The equation of 
motion is discretized on a grid of size 50x50 and solved by an adaptive-stepsize 
Runge-Kutta method. We use quasistatic driving condition, applying a field 
H and integrating the motion until the domain wall comes to rest. The system 
is first demagnetized, with the successive application of positive and negative 
fields of decreasing amplitude, and then we cycle the field between — and 
Hm- The Rayleigh parameters are extracted from the scaling properties of the 
hysteresis loops (see Fig. 1) which are averaged over several realizations of 
the disorder. The procedure is repeated for different values of the domain wall 
energy 7^,, and in Fig 2 we show that simulations are in agreement with the 
scaling theory. 



4 Random-field Ising model 

The RFIM is probably the simplest model showing the combined effect of 
disorder and exchange interaction on the shape of the hysteresis loop. The 
model can be solved exactly in one dimension and minor loops can also be 
computed [27]. Recently we have been able to obtain exact results for the 
entire demagnetization process in d — 1, including the Rayleigh laws. Here, 
we present numerical results for the Rayleigh laws in d — 2 and d — 3. 

Ind = 2, we perform a perfect demagnetization and thus obtain unambiguously 
the demagnetized state for a given realization of the disorder. This is done in 
practice changing the field by precisely the amount necessary to fiip the first 
unstable spin [24]. In this way, the field is cycled between —Hm and H„i and 
Hm is then decreased at the next cycle by precisely the amount necessary 
to have one avalanche less than in the previous cycle. This corresponds to 
decrease Hm at each cycle by an amount AH, with AH — > 0+. We thus 
obtain the demagnetizing curve and extract the Rayleigh parameters close 
to the demagnetized state for a system of linear size L = 50 (see Fig. 3) 
for different values of R. During demagnetization, in higher order loops we 
observe that the same set of spins flips at precisely the same fleld, despite the 
fact that the demagnetization proceeds and the initial state is different. 
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Performing a perfect demagnetization, however, it is not possible to go to 
very large system sizes. We have thus resorted to a different algorithm which 
perform an approximate demagnetization: instead of cycling the field between 
—Hjn and H^n we just flip the field between these two values. The no-passing 
rule [14] ensures that the system will go into the same states that would be 
obtained decreasing and increasing adiabatically H between —Hm and H^- 
The only difference is that, after each cycle, we have to decrease Hm by an 
arbitrary value AH (i.e. we can not perform the limit AH — > O"*" exactly). 
We have confirmed by numerical simulations that the results obtained with 
an approximate demagnetization with AH < 10~^ are in good agreement with 
the ones obtained under a perfect demagnetization. 

An approximate demagnetization is then used in d = 3, in order to analyze the 
effect of the phase transition on demagnetization and the Rayleigh law (the 
transition is not present in d — 1, while in d — 2 the issue is controversial). 
We find that demagnetization is possible only for R > R^, where R^. = 2.16 
for J = 1 [24]. For R < R^ the demagnetization curve coincide with the 
saturation loop and it is thus not possible to define the Rayleigh parameters. 
We thus measure a and h for different values oi R > Rc and linear system 
sizes ranging from L = 25 to L = 100. The results show that a and b vanish 
for R ^ Rc and follows a scaling law a ~ ^ RcY" and b ^ {R — Rc)^'', 
with f3a — l3b — 0.5. This result suggests that the demagnetization curve scales 
simply as M = (i? — Rc)^^"^ m{H) for if — > 0. It is interesting to compare this 
result with the behavior expected for H — He on the saturation loop: in that 
case the magnetization scales with an exponent /3 ~ 0.04, while /3 = 1/2 is 
valid only in mean-field theory. 



5 Experimental comparison and conclusions 

In order to test the result obtained above, we perform a set of experiments 
on magnctostrictive ribbons of Feg4Co2iBi5 amorphous alloy under moderate 
tensile stress [28,29]. This material is characterized by extended domain walls 
and the statistical properties of its Barkhausen noise have been recently shown 
to be well described by models of domain wall pinning [22,29]. We thus ex- 
pect that its hysteretic behavior could also be described by collective pinning 
theory. 

The sample is first demagnetized and hysteresis loops are measured at low 
fields. We find that the Rayleigh law is not perfectly verified at very low fields, 
but the hysteresis cycle becomes parabolic at higher field. At very low fields 
the magnetization is completely reversible, then for higher fields a definite and 
reproducible magnetization jump appears. Finally, at still higher fields several 
jumps combined together give rise to what can be well approximated by a 
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parabolic cycle. This behavior is summarized in Fig. 5a where we report the 
behavior of M„,jH„i as a function of i/„j. The linear behavior at high field is 
described the Rayleigh law and the parameters a and b are consistent with 
estimates obtained from the loop shape. 

In order to understand this behavior we have simulated the domain wall model 
for a single realization of the disorder using a grid of linear size L = 100. Also 
in this case, we observe at low fields a reversible response and a jump at higher 
fields. After the jump the hysteresis loop is well described by the Rayleigh law 
(see Fig. 5b). Notice that when we average the cycles over different disorder 
configurations, the jump disappears. Clearly, disorder averaging can not be 
performed directly for experimental data, unless different samples arc used. 
One should thus keep this in mind when Rayleigh parameters are estimated. 

In conclusions, we have discussed the occurrence of the Rayleigh law of hys- 
teresis, analyzing two classes of models: disordered spin models and domain 
wall depinning models. This is a step towards the ambitious goal to recover 
hysteresis loops properties from the magnetic microstructure. 

References 

[I] L. Rayleigh, Philos. Mag., Suppl. 23 (1887) 225. 

[2] G. Bertotti, Hysteresis in Magnetism, Academic Press, San Diego, 1998. 
[3] S. Chikazumi, Physics of Magnetism, Wiley, New York, 1964. 
[4] R. Vergne, J. C. Cotillard, J. L. Porteseil, Rev. Phys. Appl. 16 (1981) 449. 
[5] D. Damjanovic, J. Appl. Phys. 82 (1997) 1788. 

[6] D. Bolten, U. Bottiger, T. Schneller, M. Grossmann, O. Lohse, R. Waser, Appl. 
Phys. Lett. 77 (2000) 3830. 

[7] G. Bertotti, F. Fiorillo, M. Pasquale, J. Appl. Phys. 69 (1991) 5930. 

[8] L. Neel, Cah. Phys. 12 (1942) 1. 

[9] K. H. Pfeffer, Phys. Stat. Sol. 21 (1967) 857. 

[10] H. Kronmiiller, Z. Angew. Phys. 30 (1970) 9. 

[II] H. Kronmiiller et al., J. Magn. Magn. Mat. 13 (1979), 53. 

[12] H. Kronmiiller, T. Reininger, J. Magn. Magn. Mat. 112 (1992) 1. 

[13] A. Magni, C. Beatrice, G. Durin, G. Bertotti, J. Appl. Phys. 86 (1999) 3253. 

[14] J. P. Sethna, K. Dahmen, S. Karta, J. A. Krumhansl, J. D. Shore, Phys. Rev. 
Lett. 70 (1993) 3347. 



8 



[15] H. R. Hilzinger, H. Kronmiiller, J. Magn. Magn. Mat. 2 (1976) 11. 
[16] T. Nattermann, Y. Shapir, I. Vilfan, Phys. Rev. B 42 (1990) 8577. 
[17] Y. Enomoto, J. Magn. Magn. Mat. 129 (1994) L146. 
[18] S. T. Chui, Phys. Rev. B 51 (1995) 250. 

[19] J. S. Urbach, R. C. Madison, J. T. Markert, Phys. Rev. Lett. 75 (1995) 4694. 
[20] P. Cizeau, S. Zapperi, G. Durin, H. E. Stanley, Phys. Rev. Lett. 79 (1997) 4669. 

[21] S. Zapperi, P. Cizeau, G. Durin, H. E. Stanley, Phys. Rev. B 58 (1998) 6353. 
[22] G. Durin, S. Zapperi, Phys. Rev. Lett. 84 (2000) 4075. 

[23] S. Lemerle, J. Ferre, C. Chappert, V. Mathet, T. Giamarchi, P. L. Doussal, 
Phys. Rev. Lett. 80 (1998) 849. 

[24] O. Perkovic, K. A. Dahmen, J. P. Sethna, Phys. Rev. B 59 (1999) 6106. 

[25] L. Bocquet, H. J. Jensen, J. Phys. I (Prance) 7 (1997) 1603. 

[26] I. F. Lyuksyutov et al., Phys. Rev. B 59 (1999) 4260. 

[27] P. Shukla, Phys. Rev. E 62 (2000) 4725. 

[28] C. Appino, G. Durin, V. Basso, C. Beatrice, M. Pasquale, G. Bertotti, J. Appl. 
Phys. 85 (1999) 4412. 

[29] G. Durin, S. Zapperi, J. Appl. Phys. 85 (1999) 5196. 



9 




Fig. 1. The Rayleigh law computed in flexible domain wall in a random medium, (a) 

Hysteresis loop for different values of Hm. (b) The lower branch of the loops can be 
rescaled according to the Rayleigh law. (c) The area of the loop scales as ~ 
(the line has a slope of 3. (d) The Rayleigh parameter can be also obtained by a 
linear flt of Mm/Hm vs Hm- All units are arbitrary. 
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Fig. 3. An example of exact demagnetization process of the RFIM in d = 2 for 
R = 1.64. In the lower inset we show a detail of the main figure. In the upper inset 
we show the demagnetization curve and a small loop around the demagnetized state, 
both averaged over 30 realizations of the disorder. In the lower inset we show a detail 
of the main figure. All units are arbitrary. 
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Fig. 4. The Rayleigh parameters o and 6 are computed for different values of the 
disorder R in the RFIM in cZ = 3, the results are the average over N realization of 
the disorder: AT = 100 for L = 25. AT = 50 for L = 50, TV = 20 for L = 75 and 
A'^ = 5 for L = 100. The Rayleigh law is only found for R> and the parameters 
scale to zero as power law when R —>■ Re- All units are arbitrary. 



13 



8.5 

"§ 7.5 
H 



6.5 



5.5 



(a) 



o 



.0 



0" 



.0 







0.0005 0.001 0.0015 0.002 0.0025 

H [A/m] 



0.005 
0.004 
0.003 
0.002 
0.001 




1 ' 


1 

P " 


(b) 


P 

0' 










0' 





\ 




1 



0.5 1 1.5 2 

// [arb. units] 



Fig. 5. The peak susceptibility Mm/Hm as a function of (a) for a single realiza- 
tion of the disorder in the domain wall model and (b) for a Fe64Co2iBi5 amorphous 
ribbon. 
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